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We have studied the effect of local Coulomb interactions on the electronic structure of the molecu- 
lar magnet Mni2-acetate within the LDA+U approach. The account of the on-site repulsion results 
in a finite energy gap and an integer value of the molecule's magnetic moment, both quantities being 
in a good agreement with the experimental results. The resulting magnetic moments and charge 
states of non-equivalent manganese ions agree very well with experiments. The calculated values of 
the intramolecular exchange parameters depend on the molecule's spin configuration, differing by 
25-30% between the ferrimagnetic ground state and the completely ferromagnetic configurations. 
The values of the ground-state exchange coupling parameters are in reasonable agreement with the 
recent data on the magnetization jumps in megagauss magnetic fields. Simple estimates show that 
the obtained exchange parameters can be applied, at least qualitatively, to the description of the 
spin excitations in Mni2-acetate. 
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Introduction 

Single-molecule magnets represent a new class of 
"zero-dimensional" magnetic materials consisting of 
almost |-|nfjurL-interacting, identical point-like nanoscale 
entitiesM'clcl. These magnetic molecules have attracted a 
great deal of attention from physicists as well as chemists 
for detailed studies of many novel aspects of mesoscale 
magnetism, such as quantum spin tunneling, spin re- 
laxation of mesoscopic magnets, the implicatie 
topological spin phase (Berry phase), etc iBI 
The Mni20i2[CH3COO]i6-4H20 compoundBEl, often 
referred to as Mni2-acetate or simply Mni2, is one of 
the best-studied examples of this kind of materials. Rel- 
atively slow dynamical processes taking place in not too 
strong magnetic fields, when the intramolecular spin ex- 
citations can be ignored, have been widely studied in 
experiments, and the properties of the ground-state mul- 
tiplet S' = 10 are known very well. Correspondingly, 
most of the theoretical treatments of Mni2 have been 
based on the rigid-spin model (see e.g. Refs. p^ , p^|jr^ , 
|T6| , p7 18 19| , p0| , plp^ ) , assuming that the molecule's mag- 
netic subsystem can be represented as a single collec- 
tive spin S = 10. On the other hand, the internal mag- 
netic and electronic structure of the Mni2 molecules and 
the intramolecular interactions responsible for the for- 
mation of the molecule's collective spin have been stud- 
ied in less detail. Several experiments including inelas- 
tic neutron scatteringEj, megagauss field, magnetization 
measurementfcil, and optical absorptionE3 have been per- 
formed to clarify the internal properties of the Mni2 
molecules. The experimental data h av ji j; b een analyzed 
within the spin-Hamiltonian approachESES, however, the 
resulting information was not systematic, and often did 
not lead to unequivocal conclusions (see e.g. the discus- 



sion of different parameter sets in Ref. Pq ). 

Another approach to investigate intramolecular prop- 
erties is ah initio electronic structure calculations. So far, 
density-functional computations of the electronic struc- 
ture for Mni2 have been carried out within the local 



and generalized gradi- 
n. It is well known that for 



density approximation (LI 
ent approximation (GGA)t3l 
transition metal-oxide systems the standard density func- 
tional methods (LDA or GGA) do not give reliable infor- 
mation on tJae-enetgjv-ffaps and magnetic superexchange 
interactionsEllE3Ha'EJ, and many-body effects should be 
taken into account at least in the simplest form of the 
LDA-I-U approximationEj. This is essential, e.g., for an 
adequate description of the electronic structure and the 
spin excitation spectra in MnOED'H. Mni2-acetate is also 
an insulating metal-oxide system with an energy, 
ordered of 2 eV, whereas the existing calculation£3E3 
give a finite density of states (DOS) and small pseudogap 
at the Fermi level. As a result, the total spin per molecule 
turns out to be non-integer, in contradiction with experi- 
mental results. Futhermore, the intramolecular exchange 
couplings in Mni2 have never been calculated befoip™ 
The GGA calculations have recently been performedE£l 
for another single-molecule magnet, V15, but the result- 
ing values had to be reduced by a factor of three to 
achieve an agreement with the experiments. This dif- 
ference could also be caused by the neglect of the on-site 
Coulomb repulsion, since the latter is known to decrease 
the calculated-exchange interactions in the metal-oxide 
compoundsEOH. 

In this work, we present the results of first-principle 
LDA-I-U calculations of the electronic structure and ex- 
change interactions for Mni2 . We show that the account 
of the many-body effects within the LDA-I-U approxima- 
tion gives a reasonable picture of the electronic structure, 
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with a finite energy gap and an integer value of the to- 
tal spin of the molecule. Moreover, for the first time 
we present the calculated data for the intramolecular ex- 
change constants in Mni2. We discuss the results of our 
calculations and compare them with the existing experi- 
mental data and previous theoretical studies. 

The rest of the paper is organized as follows. In Section 
I we describe the computational scheme and the approxi- 
mations employed in this work. In Section II the resulting 
electronic structure and calculated magnetic moments in 
Mni2 are presented. Section III is devoted to the compu- 
tation of the intramolecular exchange interactions. The 
conclusions are presented in Section IV. 



I. SCHEME OF CALCULATIONS 

The calculations have heen performed for the crystal 
structure found by T. Lial2l with the simplifications sim- 
ilar to those made in Ref. ^ (all the methyl groups CII3 
have been replaced by hydrogen atoms and the water 
molecules were omitted). Therefore, the real calculations 
have been done for the system Mni20i2(IICOO)i6. Part 
of this system is sketched in Fig. |l], showing the arrange- 
ment of the manganese ions and the oxygens bridging 
them. 

To start the LDA-I-U calculations one needs to identify 
the regions where the atomic characteristics of the elec- 
tronic states have largely survived ("atomic spheres"). 
Within these atomic spheres one can expand an electron 
wave function in terms of a localized orthonormal basis 
inlma) {i denotes the site, n is the main quantum num- 
ber, I — orbital quantum number, m — magnetic number 
and a — spin index). The density matrix is defined as 



E p 



ImG 



inlrn.inlriv 



{E)dE, 



(1) 



where G"^ , . , / (E) — (inlma \ (E — H) ^1 inlm a) 
are the elements of the Green function matrix in this 
localized representation (the effective Hamiltonian H will 
be defined later). In terms of the elements of this density 
matrix the generalized LDA+U functionalL^I has 

the following form 

£;L°^+U[p^(r),{n°'}] = E^'^^^ip'^ir)] (2) 
+ E''[{n-}]-Eac[{n''}] 

Here p"^ (r) is the charge density for the electrons with the 
spin a, and E^'^^^[p''{r)] is the standard LSDA func- 
tional. Eq. (||) asserts that the LSDA suffices in the 
absence of orbital polarizations, which are given by 
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{m},cr 

-I- ((to, m" I Vee I m! ,m"') (3) 
- (to, to" I Vee I m"',TO'))n^„,n^„„,„, 



where Vee is the screened electron-electron interactions. 
The matrix elements of Vee are defined via an average 
Coulomb paranieter U and the Hund intra-atomic ex- 
change constantE3 J. Finally, the last term in Eq. (H) 
describes the correction for the double counting (in the 
absence of orbital polarizjations, Eq. (||) should reduce to 
j^LSDA^^ and is given byE3 

Edc[{n'^}] = 1-UN{N-1)-Ij[N^{N^ ^1) (4) 
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1)], 



were iV" = Tr (n^,„,) and N ^ +N-^. 

In addition to the usual LDA potential, we find an 
effective single-particle potential 



V" , ^ , 

^ mm' / J 



^ {KTO"|Ke|m',TO">;;,^,„,„ (5) 

+ ((to, to" I Vee I to', to"') 
- (to, to" I Ke |r7l'",TO')X„„,,„}- 
- c/(iv - i) + J(7V- - i) 

to be used in the effective single-particle Hamiltonian 

B = Hlsda + ^ I inlma)V^^, {inlm' a \ . (6) 

mm' 

The Linearized Muffin-Tin Orbitals (LMTO) method in 
the orthogonal approximationcj has been used for the 
LSDA calculations. 

Using our Green's function method we can obtain the 
effective exchange interaction papjneters via the energy 
variation at small spin rotationscS. The correswopding 
expression in the LDA+U scheme takes the formc2l 

/ J mm' ■^m7n'm"m"' m"m"' ' v V 

{m} 

where the spin-dependent potentials / are expressed in 
terms of the potentials of Eq. (||), 



P , = U'f , - V'^ , 

-^mm' 'mm' mm' 



(8) 



while the effective intersublattice susceptibilities (i and 
j)are defined in terms of the LDA+U eigcnfunctions ip 
as 



mm 



E "nkT "n'kj i,ilm' iJlm" i,ilm' j.jlm"" 



knn 



, EnkT — En'ki 



(9) 

As a ground state, the standard spin arrangement has 
been chosen with the spins of Mnl ions directed down 
and the spins of Mn2 and Mn3 ions directed up. This ar- 
rangement was self-consistent. To check the applicability 
of the localized spin models to Mni2 the calculations for 
the ferromagnetic configurations (all spins up) have also 
been performed. 
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II. ELECTRONIC STRUCTURE AND 
MAGNETIC MOMENTS 



Detailed calculations have been carried out for the 
three values of J7 = 4, 6, and 8 eV keeping constant 
the value of the intra-atomic Hund's exchange param- 
eter J = 0.9 eV, to check the sensitivity of the results 
with respect to variation in the local Coulomb repulsion 
U . For all values of U used in the calculations, the system 
turns out insulating with the energy gap between 1.5 — 2 
eV, and the integer moment of the ferrimagnetic ground 
state /itotai — 20 fiB- Dependence of the calculated quan- 
tities on the value of U can be inferred from the data in 
the Table |, and as one can see, for most of the proper- 
ties the dependence on U is reasonably weak. The value 
U = 8 eV seems to be slightly preferrable; it is known 
that this value gives good agreement with experiments 
for other manganese-oxide systems, and agrees reasoiu 
ably well with the independent theoretical calculationEZi 
giving the value oi U = 6.9 eV. Therefore, below we dis- 
cuss in detail only the case of J7 = 8 eV. 

To begin with, let us compare the results given by LDA 
and LDA-I-U approaches. In Figs. ^ and ||, we show the 
calculated partial densities of states (DOS) for different 
manganese atoms and for the oxygen atoms belonging to 
the exchange bridges. Our LDA curves for the DOS are 
similar to those obtained in Refs. |2^,|2^ using the GGA 
approach. The main distinction of the LDA-I-U results 
as compared with the LDA is the nonzero energy gap 
instead of the finite DOS at the Fermi level. Moreover, 
there is a difference in the DOS peaks structure. All 
manganese ions in Mni2 molecule are in distorted oc- 
tahedral coordination, correspondingly, one can see the 
three peaks of t2g character situated below the two Cg- 
like peaks. For the LDA calculations the spectral density 
is strongly concentrated within the region ±2 eV near 
the Fermi level, see Fig. ^, while for the LDA-I-U calcula- 
tions, the DOS is spread rather uniformly over the region 
of about 10 eV. 

It is known that the LDA-I-U approach is effective for 
consideratioB of charge separation and charge ordering 
phenomenaEJ. Mni2 molecules also exhibit Mn'^+ - Mn^+ 
disproportionality, which has been directly demonstrated 
by the measurements of the X-ray absorption spectraE^. 
The majority spin DOS curves (see Fig. g) demonstrate 
three occupied t2g peaks and two empty eg peaks for 
Mnl, while for Mn2 and Mn3 ions only one Cg peak is 
empty; this picture indeed corresponds to the charge dis- 
proportionality detected in the experiments. 

Several remarks are in order. There are indications 
that the manganese-oxygen bonds are considerably cova- 
Icnt, and 3d states of manganese are highly hybridized 
with 2p states of oxygen. As can be seen in Fig. |[ the 
main features in the partial DOS of manganese ions have 
counterparts in the partial DOS of oxygen ions. Another 
indication of the strong covalence is the discrepancy be- 
tween the sum of the local 3c?-moments of Mn ions and 
the total magnetic moment of the molecule which is cal- 



culated from the total DOS. The difference, which is of 
order of 2fj,B, can be attributed to the partial moments 
of oxygens estimated to be of order of O.OS^b. 

The value of the energy gap calculated within LDA-j-U 
scheme varies from 1.78 eV to 2.01 eV with the change 
in the value q£ U from 6 eV to 8 eV. Recent optical 
measurementaHj show the features in the absorption and 
reflection spectra in the region of the charge-transfer ex- 
citations at 1.2-|eV, 2.0 eV, and 2.4 eV. Magneto-optical 
measurementalj also show the latter feature. These val- 
ues are in very good agreement with the calculated en- 
ergy gap. 



III. INTRAMOLECULAR EXCHANGE 
INTERACTIONS 

At present, there is no reliable quantitative informa- 
tion about the intramolecular exchange couplings in the 
Mni2 molecules, making difficult the validation of the 
corresponding computational results. Crude estimates 
for the exchange couplings, based on the comparison 
with different manganese-oxygen compounds, have been 
presented in Ref. Quantitative evaluation of the ex- 
change constants has been performed in Ref. 2^ using the 
results of the megagauss-field experiments, however the 
magnetization curves appeared to be not very sensitive to 
the values of the exchange parameters. Another quanti- 
tative approach, based on the inelastic neutron scattering 
data, has been undertakencJ, but a simplified eight-spin 
model was used which makes the fitting not conclusive 
(see the discussion of various possible sets of parameters 
in Ref. ^6|). Therefore, we believe that it makes sense to 
focus the discussion of our results on qualitative trends in 
the calculated parameters instead of their exact values. 

The exchange parameters have been calculated for dif- 
ferent values of U, see Table |. As expected for the su- 
perexchange mechanism, the values of the couplings de- 
crease with increasing U; again, below we discuss only 
the values corresponding to C/ = 8 eV. We have com- 
pared the calculated values of the exchange parameters 
with the data suggested in Refs. see the Table 
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In this Table, we have shown only the interactions 
between relatively closely situated manganese ions; the 
exchange parameters between distant ions do not exceed 
2-3 K. As one can see, the values obtained in our cal- 
culations are closer to those resulting from the analysis 
of the megagauss-field experiments, while the data sug- 
gested from general chemical consideration are bigger by 
about a factor of three. However, the more interesting 
fact is that the exchanges between the Mn^'^ ions are 
not all equal, as assumed in Refs. The Mn'*+ ions 

in Mni2 molecules are arranged in the four corners of a 
cubane (distorted cube, see Fig. ||) while the other four 
corners are occupied by the oxygens providing superex- 
change bridges between the manganese ions. Therefore, 
the manganeses form a distorted tetrahedron with four 
short edges of the length 5.33 a.u. and two long edges of 
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the length of 5.56 a.u., and the angles Mn-O-Mn in the 
cubane which correspond to the long and short edges of 
the manganese tetrahedron are 95.5° and 100.6°. In spite 
of the moderate difference in the angles and lengths of 
the exchange pathways, the calculated coupling param- 
eters turn out very different; the interactions J3 for the 
longer pathways are about three times smaller than the 
interactions J3 for the shorter ones_-||— . 

All the previous considerationscjEScJ of the exchange 
interactions in Mni2 were based on the localized Heisen- 
berg model, assuming that the values of the magnetic 
moments and the exchange parameters do not depend 
on the spin arrangement of the molecule. To check the 
accuracy of this assumption, we have performed the cal- 
culations for the totally ferromagnetic spin configuration, 
when all manganese spins are directed up, using the pa- 
rameters U = 8 eV and J — 0.9 eV. The results are the 
following: 

/^Mnl — 3.06 flB, AiMn2 = 3.55 ^B, ^J-Mn3 = 3.85 /i_B , 

gap = 1.63 cV, 

Ji = -36 K, J2 = -18 K, 

J3 -22 K, J3 = -1-1 K, Ji = -6 K. 

Comparing these results with the Tables |l III 



we see 



that the values of the moments vary by no more than 5%, 
but the exchange integrals Ji and J3 decrease by about 
25-30%, and the exchange J3 almost vanishes. These 
changes are due to the strong covalency effects and the 
corresponding partial delocalization of the 3c?-electrons 
discussed above. 

The energy of the exchange excitations in the Mni2 
molecule, e.g., the difference between the ground-state 

= 10 multiplet and the close excited multiplet, can be 
crudely obtained from a simple mean-field estimate, as- 
suming that the pairs Mn'*"*'-Mn'^"'' ions form stiff dimers 
with the total spin 1/2 (i.e., Ji —00). Then, the en- 
ergy of the ground state is Eq — — 2 J3 — J3 — 32 J4 + 16 J2 
and the energy of the first excited state is Ei = 2 J3 -f J3, 
so that the excitation energy is 52 K. This number is close 
to the values-58-66 K obtained in the neutron scattering 
experimentaHj. Of course, this agreement should not be 
taken literally (it is a semiclassical estimate, the value 
of Ji is not big enough, we did not take into account 
the contributions of relativistic interactionS|-ss[hich could 
be comparable to the exchange contributioncB, etc.), but 
we believe that this indicates qualitative adequacy of the 
calculated exchanges. 



IV. CONCLUSIONS 

We have applied the LDA-I-U approach to the study 
of intramolecular magnetic interactions in Mni2-acetate 
system. We have shown that the account of the on-site 
Coulomb interactions provides a good description of the 
main features of the electronic and magnetic structure, 
giving the finite energy gap and the integer value for 
the magnetic moment of the molecule. The calculated 
partial magnetic moments, charge disproportionality be- 
tween the Mnl and Mn2, Mn3 ions, and the energy gap 
are in agreement with existing experimental data. The 
calculations show significant p-d hybridization and con- 
siderable covalency effects, and the contribution of de- 
localized electrons in the total magnetic moment of the 
molecule is about 10%. The calculated intramolecular ex- 
change interactions are in qualitative agreement with the 
existing data. An interesting result of the calculations is 
that the interactions between the Mn^+ ions situated in 
the inner cubane are not all equal, differing by about a 
factor of three for different exchange-coupled pairs. The 
applicability of localized Heisenberg model has been esti- 
mated by comparing the values of the exchange interac- 
tions for different spin configurations (ground-state ferri- 
magnetic and totally ferromagnetic), and corresponding 
variation of the coupling parameters is about 30%. 
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FIG. 1: The sketch of the arrangement of the manganese and 
oxygen ions in the Mni2 molecule. Big dark-grey circles rep- 
resent the Mn'*"'" ions, big light-grey circles — Mn^"*" ions, and 
small dark circles are the bridging oxygens. Non-equivalent 
manganese ions Mnl, Mn2, Mn3, and the oxygens 01, 02, 
03 bridging them are indicated. 
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FIG. 2: Partial densities of states for different manganese (d 
orbitals) and oxygen {p orbitals) atoms within the LDA+U 
approach. The values Z7 = 8 eV and J = 0.9 eV have been 
used. 



TABLE I: The calculated magnetic moments (in /is) of the 
manganese ions, values of the exchange interaction parame- 
ters (in K) and the energy gap (in eV) for different values of 
U (in eV). The value of J is 0.9 eV. The parameters of the 
exchange interactions J3 correspond to the exchange between 
the closest Mnl ions in the cubane (see text). 
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FIG. 3: Partial densities of states for different manganese (d 
orbitals) and oxygen (p orbitals) atoms calculated within the 
LDA approach. 



TABLE II: The comparison of the calculated magnetic mo- 
ments (in /is) with the existing theoretical and experimental 
data. 1 are the data calculated by Poderson et al., Ref. [26,27]; 
2 — the data calculated by Zeng et al., Ref. [28] (only the 
contribution of d-states is presented) ; Exp — the experimen- 
tal results of Robinson et al., Ref. [41]. The next two columns 
are the results of our calculations done within the LDA and 
LDA-I-U methods. LDA-I-U results for the values U = 8 eV, 
J = 0.9 eV are shown. 
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FIG. 4: The sketch of the arrangement of the manganese and 
oxygen ions in the inner cubane of Mni2 molecule. Arrows 
show the directions of the displacement of the Mn ions which 
lead to the formation of the inner cubane instead of the perfect 
cube. 



TABLE III: The comparison of the calculated values of the 
exchange parameters (last line of the table) with the existing 
data. 3 — the data suggested by Scssoli ot al., Ref. [43]; 4 
— the data of Platonov ct al., Rcf. [23]. LDA+U results for 
the values U = 8 eV, J = 0.9 eV are shown. The parameters 
of the exchange interactions Js and J3 correspond to the ex- 
change inside different pairs of Mnl ions in the cubane (see 
text). All values are given in Kelvin. 
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